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The capability to establish adaptive relationships with the environment is an essential 
characteristic of living cells. Both bacterial computing and bacterial intelligence are 
two general traits manifested along adaptive behaviors that respond to surrounding 
environmental conditions. These two traits have generated a variety of theoretical and 
applied approaches. Since the different systems of bacterial signaling and the different 
ways of genetic change are better known and more carefully explored, the whole adaptive 
possibilities of bacteria may be studied under new angles. For instance, there appear 
instances of molecular "learning" along the mechanisms of evolution. More in concrete, 
and looking specifically at the time dimension, the bacterial mechanisms of learning 
and evolution appear as two different and related mechanisms for adaptation to the 
environment; in somatic time the former and in evolutionary time the latter. In the 
present chapter it will be reviewed the possible application of both kinds of mechanisms 
to prokaryotic molecular computing schemes as well as to the solution of real world 
problems. 
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INTRODUCTION 

Bacterial computing, as an applied field recently launched (Poet 
et al., 2010), as well as the theoretical approaches to prokaryotic 
or bacterial intelligence, are derived from the adaptive response of 
living cells to existing environmental conditions. From a practical 
standpoint, we could define bacterial computing as the possibil- 
ity of using bacteria for solving problems that today are solved 
by computers. If a bacterium could perform the work of a com- 
puter, this would allow us to build millions of computers which 
be replicated every 30min, and that they would be confined 
within a Petri dish. According to Amos (201 1) natural computing 
paradigms inspired by biological processes (e.g., artificial neural 
networks, genetic algorithms, ant colony algorithms, etc.) have 
proved to be very effective. However, all these "forms of comput- 
ing" occurs in silico, and therefore within a computer. At present 
and in agreement with Amos (2011), the challenge is the pos- 
sibility to use biological substrates and biological processes to 
encode, store and manipulate information (Cordero et al, 2013). 
For instance, to build a simple computing device, using bac- 
teria rather than silicon. Since the seminal work of (Adleman, 
1994) the feasibility of using biological substrates for computing 
has been well-established: Levskaya et al. (2005) has shown that 
the living cell could be considered as a programmable compu- 
tational device, Baumgardner et al. (2009) using DNA segments 
and Hin/hixC recombination system successfully programmed E. 
coli with a genetic circuit that enables bacteria to solve a classi- 
cal problem in artificial intelligence, the Hamiltonian problem; or 
the theoretical model where bacteria are used to solve the "burnt 
pancake problem" (Heyer et al, 2010). 

In a theoretical realm, bacterial computing could be an 
emergent phenomenon consequence of learning and evolution. 



Bacterial learning and evolution are but two different and related 
mechanisms for adaptation to the environment, in somatic time 
the former and in evolutionary time the latter (Di Paola et al., 
2004). In this chapter we review the possible application of both 
mechanisms to prokaryotic molecular computing as well as to the 
solution of real world problems. 

During recent years, some experiments have shown that 
bacteria can learn the ability to anticipate changes in their 
immediate environment. For instance, Tagkopoulos et al. (2008) 
found how E. coli colonies can develop the ability to associate 
higher temperatures with a lack of oxygen, and how bacteria 
have naturally "learned" to get ready for a serving of maltose 
after a lactose appetizer (Mitchell et al., 2009). According to 
Tagkopoulos et al. (2008), homeostasis explains microbial 
responses to environmental stimuli — by means of intracellular 
networks, microbes could exhibit predictive behavior in a fashion 
similar to metazoan nervous systems. Even more, bacteria are 
able to explore the environment within which they grow by 
utilizing the motility of their flagellar system (Lahoz-Beltra, 
2007) and deploying a sophisticated "chemotactic" navigation 
system that samples the environmental conditions surrounding 
the cell and systematically guides away from the unfavorable 
conditions and toward the favorable ones. 

In this chapter we review several theoretical studies, mod- 
els, and simulations about bacterial forms of natural computing, 
gauging its potential application and impact. We review how 
proteins form molecular complexes and networks related to 
molecular signaling functions and bacterial information pro- 
cessing. Modeling proteins as McCulloch- Pitts neurons reviews 
a hardware model (Di Paola et al, 2004) demonstrating how 
proteins of the intervening signal transduction networks could be 
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modeled as artificial neurons, simulating the dynamical aspects 
of the bacterial taxis. The model is based on the assumption that, 
in some important aspects, proteins (Di Paola et al, 2004) of the 
signaling system may be considered as McCulloch-Pitts artificial 
neurons (McCulloch and Pitts, 1943) that transfer and process 
information from the bacterium's membrane to the flagella motor. 

Modeling proteins as networks of processing elements will 
review how these proteins are also involved in other bacterial 
signaling functions through complex molecular systems. Finally, 
Modeling metabolites as "metabolic hardware" suggests how sim- 
ilar "informational" properties of proteins, particularly enzymes, 
may organize cellular metabolism, by introducing in our study 
the concept of "metabolic hardware." 

Bacteria can also evolve some true learning behaviors to 
respond optimally to their environment (Bacterial evolution). 
At present, most methods in evolutionary computation are 
inspired in the fundamental principles of neo-Darwinism or 
population genetics theory, considering as main sources of 
variability chromosome crossover (or recombination) and muta- 
tion (Per ales -Gravan and Lahoz-Beltra, 2008). However, bacteria 
exhibit several other genetic mechanisms as sources of variabil- 
ity, i.e., mechanisms such as transformation, conjugation, and 
transduction. In Bacterial conjugation we discuss the applicabil- 
ity of conjugation, a genetic mechanism exhibited by bacterial 
populations, and we simulate the evolutionary process along 
this mechanism. The efficiency of bacterial conjugation (Perales- 
Gravan and Lahoz-Beltra, 2008) is illustrated designing, by means 
of a genetic algorithm based on this very mechanism, an AM radio 
receiver. In Bacterial transduction we continue the study of bacte- 
rial evolution, in this case by modeling and simulating the whole 
bacterial transduction mechanism. 

MODELING PROTEINS AS McCULLOCH-PITTS NEURONS 

Adaptive behavior in bacteria depends on organized networks of 
proteins (Figure 1) governing molecular processes within the cel- 
lular system. Bacteria are able to explore the environment within 
which they develop by utilizing the motility of their flagellar 
system as well as a biochemical navigation system that sam- 
ples the environmental conditions surrounding the cell (Di Paola 
et al., 2004). In this chapter we described how in some impor- 
tant aspects proteins can be considered as processing elements 
or McCulloch- Pitts artificial neurons that transfer and process 
information from the bacterium's membrane surface to the flag- 
ellar motor. The McCulloch- Pitts artificial neuron (Lahoz-Beltra, 
2008) is a mathematical model (Di Paola et al, 2004) of a bio- 
logical neuron. A neuron has a set of inputs JiJ^v . . Jb a set of 
weight values associated with each input line Wi,W2,. . . ,Wj, one 
output Of and a linear threshold function f(netj). Every time a 
neuron receives a set of input signals, performs the weighted sum 
(with the weights associated with each line) obtaining a neti value, 
finally deciding its state or output with a threshold function: 



CW CCW 



neti 



1 neti > 1 
0 neti < 1 



(i) 



(2) 




■> CheY ) 
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FIGURE 1 | Che protein network mediating bacterial taxis. External 
factors are identified by membrane receptors (MCP) which reactivity 
depends on methylation levels controlled by CheB and CheR cytoplasmic 
proteins. In the network CheW is the bridge between the receptors and 
protein kinase CheA which can donates a phosphate group to CheY. The 
binding of CheY-P to the switch complex protein FliM induces a clockwise 
flagellar rotation (CW) whereas the dephosphorylation of CheY (via CheZ) 
restores counter-clockwise flagellar rotation (CCW) (Di Paola et al., 2004). 



The hardware model of the McCulloch-Pitts output artificial 
neuron assumes some plausible analogies between neurons and 
proteins. For instance, the connection between neurons, thus 
synapses as well as the neuron input and output would have 
their equivalent in the proteins on the concepts of bond strength, 
external factors (e.g., heat, ions, chemical substances, etc.) and 
conformational states related to catalysis and binding, respec- 
tively. In this analogy, it is also assumed that activation function 
represents cooperativity in proteins. 

The hardware model of the proteins as a McCulloch- Pitts 
artificial neuron is based on a 741 inverting operational ampli- 
fier (Rietman, 1988). The operational amplifier simulates the 
bond strength, external factors, conformational states and coop- 
erativity by means of a potentiometer (Ri,R2,...,R n ), voltage 
(Vi,V2,...,V n ), voltage (Vo) and its saturation curve. In the cir- 
cuit Ri,R2,...,R n are n variable resistors modeling the weights 
associated to connections among the n membrane receptors — 
playing the role of the input neurons — and the McCulloch- Pitts 
output neuron. Thus, variable resistors simulate the degree of 
influence of n inputs or external factors acting as repellents assum- 
ing one input per membrane receptor. The resistor Rf simulates 
the feedback in the output neuron being Vo the output voltage. 
Assuming the presence in the medium or environment of n inputs 
or external factors which presence is modeled as Vi(i = 1, . . . , n) 
input voltages with Ri(i = 1, . . . , n) weights or variable resistors. 
According this model (Di Paola et al, 2004) it follows that: 



where 0 is a threshold value. 



V 0 = -R f J2 



i=l 



Vi 

Ri 



(3) 
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FIGURE 2 | (A) Hardware and (B) software operational amplifier circuits 
with three repellent inputs. A voltage of +3 V is applied to the 
potentiometers while the amplifier is powered with +9V (+4.5 and —4.5 V). 
In the hardware circuit (A) a LED diode connected to the amplifier output 
detects the intensity of the device signal. In the software circuit (B) a virtual 
voltage probe replaces to the LED diode (Di Paola et al., 2004). 



Once the input is applied, the value of the output is given by the 
LED diode state powered by the Vb voltage. Even though LED 
diode intensity changes according to the output voltage in our 
experiments we only considered two states of the LED diode. 
When LED diode state is switch- off state then it means an out- 
put equal to 0 simulating that CheY and FliM remain separated 
and by consequence the flagellar rotation is CCW, indicating the 
bacterium swimming behavior. Otherwise, when the LED diode 
is switched-on then it means an output equal to 1, simulating 
that CheY binds to FliM, and as a result the flagellar rotation 
is CW, indicating a bacterium tumbling behavior. The thresh- 
old between states simulates the critical level of phosphorylation 
in which CheY binds to FliM resulting in the transition from 0 
(swimming behavior) to 1 (tumbling behavior). 

In the practical implementation of the model, the amplifica- 
tion factor A of the operational amplifier: 

A =- R fibj (4) 

i=i j 

simulates the catalytic effect of enzymes, e.g., the protein kinase 
CheA which donates a phosphate group to CheY resulting a phos- 
phorylated CheY. As a consequence, if we establish a similarity 
between the Rf resistor and the Michaelis-Menten K m constant 
then the resistor Rf would be simulating the affinity between 
the enzyme and substrate. Note that in our hardware model 
(Figure 2) the inhibitory connection that is usually included 
in the hardware implementation of artificial neural networks 
(Rietman, 1988) has been removed. 

MODELING PROTEINS AS NETWORKS OF PROCESSING ELEMENTS 

In general, the molecular systems involved in bacterial signal- 
ing (and in M. tuberculosis) are extremely diverse, ranging from 
very simple transcription regulators (single proteins compris- 
ing just two domains) to the multi- component, multi-pathway 
signaling cascades that regulate crucial stages of the cell cycle, 
such as sporulation, biofilm formation, dormancy, pathogene- 
sis, etc. (Navarro and Marijuan, 2011). A basic taxonomy of 
bacterial signaling is shown in Figure 3. The first level of complex- 
ity corresponds to the simplest regulators, the "one-component 
systems." Actually, most cellular proteins that participate in cel- 
lular adaptation to the changing environment, in a general sense, 
could be included as participating within this elementary cate- 
gory (Galperin et al, 2001). Following the complexity scale is 
the "two-component systems," which include histidine protein- 
kinase receptors and an independent response regulator; they 
have been considered as the central signaling paradigm of the 
prokaryotic organisms, since a number of intercellular and inter- 
species communication processes are served by these systems. A 
further category (conceptual consistency) of "three-component 
systems" is applied to those two -component systems that incor- 
porate an extra non-kinase receptor to activate the protein-kinase 
(Marijuan et al, 2010). 

The Signaling/Transcriptional Regulatory Networks, like the 
M. tuberculosis network made for the authors of this paper (Sanz 
et al., 201 1), may be used to analyze the ability of Mycobacterium 
to perceive the host signals in different tissues and cell types, as 
well as adaptive responses that bacteria organizes against them. 



In that sense, to adequately study on processes of latency and 
reactivation using these networks would be very important. These 
networks will be useful to provide an overview of multiple func- 
tional aspects of this bacterium and to suggest new experiments. 

MODELING METABOLITES AS "METABOLIC HARDWARE" 

In this section we review the possibility of using metabolic net- 
works as hardware in the study of the optimization of metabolic 
pathways as well as in the field of molecular and natural com- 
puting. We call to these bioinspired architectures as metabolic 
hardware. In particular, adopting as an example well known 
metabolic pathway of Krebs cycle, we introduce the methodology 
(Recio Rincon et al, 2013) to translate the molecular struc- 
ture or topology of their metabolic intermediates to a binary 
matrix, showing how sugars and other glycolytic molecules could 
be modeled as binary matrices as well as LED dot matrices. In 
prokaryotic cells and bacteria which lack mitochondria, the Krebs 
cycle is performed in the cytosol. 

From a historical perspective one of the first procedures to 
translate the molecular topology to a matrix was introduced by 
Randic (1974), taking an element a-^ the value 1 when the vertices 
are adjacent or 0 otherwise. Figure 4 illustrates an example of 
this method for vitamin A or retinol (Lahoz-Beltra, 2012a). Our 
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FIGURE 3 | The one-two-three Component Systems. These three systems 
are the characteristic classes of signaling pathways developed by prokaryotes. 
The external stimulus is perceived either by an internal receptor-transducer 
(left), or by a transmembrane histidine kinase that connects with a response 



regulator (center), or by an independent receptor associated to the histidine 
kinase (right). This scheme represents the basic taxonomy of bacterial signaling; 
the three different options imply very different information processing 
capabilities and metabolic costs (Modified from Marijuan et al., 2010). 



method assigns a 5 -bit word to the functional groups of the 
molecule (Figure 5). For that purpose we define a table or Rosetta 
stone (Table 1) that includes the most frequent functional groups 
in metabolic intermediates, which were ordered by its redox 
potential (tendency of a functional group to acquire electrons). 

Let S and P be two binary matrices which represent respec- 
tively the substrate S m and product P m of a biochemical reaction 
catalyzed by an enzyme E m . Since in Krebs cycle all metabolites or 
metabolic intermediates are molecules with 4 or 5 carbon atoms, 
we defined (5) and (6) matrices respectively: 



C 4 : 



c 5 = 



( an an an an «is \ 

#21 a 22 a 23 &24 <?25 
#31 #32 #33 #34 #35 
\ 041 #42 043 044 045 / 

/ an an an an an \ 

#21 #22 #23 #24 #25 

#31 #32 #33 #34 #35 

041 042 043 #44 #45 

\#51 #52 #53 #54 #55 / 



(5) 



(6) 



Note that given a value i, (an, an, ... , 0/5) is a row vector 
representing the functional group of the substrate Sy or product 
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FIGURE 4 | The molecule of vitamin A or retinol represented as a 


binary matrix (Lahoz-Beltra, 2012a; Transl: Polish). 
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FIGURE 5 | Glucose molecule (Left) and its hardware version as a matrix of LEDs (Right) simulated with CEDAR Logic Simulator (Recio Rincon et al., 
2013). 



pij molecules. Thus, each row in the matrices C4 and C5 repre- 
sents a carbon atom in the molecule, having a total of 32 possible 
binary vectors from 00000 to 11111 (Table 1). Using as a criterion 
the redox potential vectors were classified from its most reduced 
(addition of hydrogen or the removal of oxygen) form or alkyl 
group to the most oxidized (addition of oxygen or the removal of 
hydrogen) or CO2. However, since the metabolites of Krebs cycle 
are the result of assembling functional groups among a total of 
22 combinations of carbon, then 10 binary vectors are without 
chemical meaning. In order to perform future simulation exper- 
iments, molecules of CO2 and acetyl- Co A were represented as a 
row vector (7) and 2x5 matrix (8) shown below: 



C0 2 = 
acetyl- Co A = 



(mi 1), 

11101 
01000 



(7) 
(8) 



Using this method the Krebs cycle was modeled as follows 
(Recio Rincon et al, 2013): 
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where each matrix stands for one of the following metabolites: 
Citrates Iso-citrate —> a-Ketoglutarate^ Succinyl-CoA —> 
Succinate —> Fumarate —> Malate —> Oxalacetate. 

BACTERIAL EVOLUTION 

At present, all methods in Evolutionary Computation (genetic 
algorithms, evolutive algorithms, genetic programming, etc.) are 
bioinspired by the fundamental principles of neo-Darwinism 
(Lahoz-Beltra, 2008), and by a vertical gene transfer; that is to 
say, by a mechanism in which an organism receives genetic mate- 
rial from the ancestor from which it evolved (Per ales -Gravan 
et al, 2009). Indeed, most thinking in Evolutionary Computation 
focuses upon vertical gene transfer as well as upon crossover 
and/or mutation operations. 

Bacteria are microscopic organisms whose single cells repro- 
duce by means of a process of binary fission or of asexual 
reproduction, bearing a resemblance to John von Neumanns 
universal constructor (Von Neumann, 1966). Thus, a bacte- 
rial population (or colony) evolves according to an evolutive 
algorithm similar to Dawkin's biomorphs (Dawkins, 1986), the 
cumulative selection of mutations powering their evolution. 
Bacteria, however, exhibit significant phenomena of genetic 
transfer and crossover between cells. This kind of mechanism 
belongs to a particular kind of genetic transfer known as hor- 
izontal gene transfer. Horizontal, lateral or cross-population 
gene transfer is any process in which an organism, i.e., a donor 
bacterium, transfers a genetic segment to another one, a recipient 
bacterium, which is not its offspring. In the realm of biology, 
whereas the scope of vertical gene transfer is the population, 
in horizontal gene transfer the scope is the biosphere. This 
particular mode of parasexuality between "relative bacteria" 
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includes three genetic mechanisms: conjugation, transduction 
and transformation. Furthermore, microorganisms are very 
interesting individuals because they also exhibit "social interac- 
tions." We found (Lahoz-Beltra et al., 2009) how the inclusion 
of the "social life of microorganisms" into the genetic algorithm 
cycle, significantly improves the algorithm's performance. 

In this section we explore the possibility of using for practical 
purposes some of the observed genetic transfer mechanisms in 
bacteria. 

BACTERIAL CONJUGATION 

This section describes a biologically inspired conjugation opera- 
tor simulating a bacterial conjugation. Its usefulness is illustrated 
in a set of computer simulation experiments where including 
such operator into a genetic algorithm we were able to design an 
AM radio receiver (Perales-Gravan and Lahoz-Beltra, 2008). The 
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FIGURE 6 I Homologous recombination or cross over mechanisms. (A) 

One-point re-combination. (B) Two-points recombination. Bacterial 
conjugation and recombination. (C) With a random point (CORP). (D) With a 
fixed point (COFP) on the donor chromosome (Perales-Gravan and 
Lahoz-Beltra, 2008). 



attributes optimized by this algorithm include the main features 
of the electronic components of an AM radio circuit, as well as 
those of the radio enclosure designed to house the radio circuit 
(Perales-Gravan and Lahoz-Beltra, 2008). 

A bacterial genetic algorithm is an evolutionary strategy based 
on bacterial conjugation and mutation. Starting with a random 
population of circular chromosomes reproduction, conjugation 
and mutation were simulated, obtaining new generations of equal 
size. The current bacterial algorithm uses homologous recom- 
bination after conjugation, a population size and a conjugation 
parameter, as well as a conjugation (or recombination) and muta- 
tion probabilities (Perales-Gravan and Lahoz-Beltra, 2008). Note 
that in the biological realm as well as in the simulations, the term 
population could be substituted by strain or colony and the lin- 
ear chromosomes of a genetic algorithm are replaced by circular 
chromosomes. Our operator which includes the recombination 
between bacterial chromosomes assumes that donor bacterium is 
always Hfr. Two different conjugation operator versions (Perales- 
Gravan and Lahoz-Beltra, 2008) have been defined (Figure 6). In 
both definitions since transfer of the donor bacterial chromosome 
is almost never complete, then the length of the strand transferred 
to the recipient cell has been simulated applying Monte Carlo's 
method and assuming DNA lengths exponentially distributed: 



/ = ln(U) 

a 



(9) 



being U a random number and a the conjugation parameter. 
The conjugation parameter summarizes all the relevant factors 
affecting the length / value. One of the most relevant factors 
affecting value is the temperature promoting the agitation of the 
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bacteria, disrupting conjugation before the entire chromosome 
can be transferred (Perales-Gravan and Lahoz-Beltra, 2008). 

In bacteria, crossing over involves the aligning of the donor 
chromosome segment with its homologous segment on the 
recipient bacterial chromosome. Next, a break occurs at a point 
origin and an end point of the recipient chromosome, removing 
and replacing the segment with corresponding homologous genes 
from the segment of the donor chromosome (Perales-Gravan 
and Lahoz-Beltra, 2008). The described steps are repeated several 
times, thus a number of times equal to the bacteria population 
size. The efficiency of the bacterial conjugation operator has been 
illustrated designing an AM radio receiver with a genetic algo- 
rithm based on this operator (Perales-Gravan and Lahoz-Beltra, 
2008). 

In Figure 7, we show the bacterial chromosome coding for 
the main features of the radio receiver and in Figure 8 a 



representative performance graph (average fitness per generation) 
of the experiments where simulated bacterial colonies evolved 
searching for the optimized circuit and enclosure. 

BACTERIAL TRANSDUCTION 

In Nature, microorganisms such as bacteria and viruses share 
a long and common evolutionary relationship (Perales-Gravan 
et al., 2009). This relationship is mainly promoted by bacterio- 
phages (or phages), a kind of virus that multiplies inside bacteria 
by making use of the bacterial biosynthetic machinery. Some 
bacteriophages are capable of moving bacterial DNA (the "bacte- 
rial chromosome") from one bacterium to another. This process 
is known as transduction. When bacteriophages infect a bacte- 
rial cell, their normal mode of reproduction makes use of the 
bacterium's replication machinery, making numerous copies of 
its own viral genetic material (i.e., DNA or RNA). The nucleic 
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FIGURE 7 | Bacterial chromosome (left) with 14 genes codifying for the main characteristics of an AM radio receiver (right) (Perales-Gravan and 
Lahoz-Beltra, 2008). 
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FIGURE 8 | Performance graph where bacterial colonies evolved searching for an optimized AM radio receiver. (A) CORR (B) COFP (Perales-Gravan and 
Lahoz-Beltra, 2008). 
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acid copies (or chromosome segments) are then promptly pack- 
aged into newly synthesized copies of bacteriophage virions. 
Generalized transduction occurs when "any part" of the bacte- 
rial chromosome (rather than viral DNA) hitchhikes into the 
virus (i.e., T4 phages in E. coli bacterium). However, when only 
"specific genes" or certain special "segments" of the bacterial 
chromosome can be transduced, such a mistake is known as 
specialized transduction (i.e., X phages in E. coli bacterium). 

In transduction, transference of chromosome segments 
between bacterial populations or colonies is very different from 
migration (the occasional exchange of individuals). Migration 
and transduction could bear a resemblance, but only when trans- 
duction involves the complete chromosome transference between 
bacterial populations. Furthermore, this kind of transference is 
a highly unlikely event in bacteria, transduction of chromosome 
segments taking place in these microorganisms. 

In this section, we model and simulate the two kinds of 
transduction operations (Figure 9) examining the possible role 
and usefulness of this genetic mechanism in genetic algorithms 
(Perales-Gravan et al., 2013). In a previous section (Perales- 
Gravan and Lahoz-Beltra, 2008), we introduced a bacterial 
conjugation operator showing its utility by designing an AM 
radio receiver. Conjugation is one of the key genetic mech- 
anisms of horizontal gene transfer between bacteria. In the 
present section, we refer to a genetic algorithm including trans- 
duction as PETRI (Promoting Evolution Through Reiterated 
Infection). We investigated the transfer of genes and chromo- 
somes among sub-populations with a simulated "bacteriophage." 
In the model we consider a structured population divided among 



several sub-populations or "bacterial colonies," bearing a resem- 
blance with coarse-grain distributed genetic algorithms. Each 
sub-population is represented as a Petri dish (a glass or plastic 
cylindrical dish used to culture microorganisms). It should be 
noted, however, that even when we divide a population into sub- 
populations, the proposed algorithm is sequential. Thus, the algo- 
rithm is not a distributed one, since we used a mono-processor 
computer and the algorithm was not parallelized. Moreover, the 
migration mechanism is synchronous, as gene and chromosome 
transferences were both between sub-populations and during the 
same generation. Therefore, our approach could be related with 
those models of Cellular Genetic Algorithms (cGA) adopted also 
for mono-processor machines (Alba and Dorronsoro, 2008), with 
no relation to parallelism at all. In our model, we assumed that 
bacteria are capable of displaying crossover through conjugation, 
instead of performing one-point or two-point recombination. 
Moreover, we assume that no vertical gene transfer mechanism 
is present in bacterial populations. 

With the aim of studying the performance of the transduction 
operator, we used different optimization problems. Experiments 
conducted in the presence of transduction were compared with 
control experiments, performed in the absence of transduction. 
Similarly, we compared the transduction performance under the 
three types of crossover: conjugation, one-point or two-point 
recombination. We are interested in the study of genetic algo- 
rithms based on horizontal gene transfer mechanisms, mainly 
conjugation and transduction operations. It is important to note 
that even when conjugation and transduction are both horizon- 
tal gene transfer mechanisms, there are some relevant differences 
between both. In the first place, whereas conjugation involves two 
bacteria from the same population, the bacteria involved in trans- 
duction can belong to different populations. As a consequence, 
conjugation is a genetic mechanism of horizontal gene transfer 
within a population, whereas transduction is a genetic mecha- 
nism of horizontal gene transfer between populations. Secondly, 
in conjugation, the length of the transferred genetic segment is 
variable, whereas in transduction, the transferred segment length 
is always constant. 

Let b be a chromosome (i.e., bacterium; 1, . . . , j, . . . , N) and 
p a sub-population (i.e., Petri dish; 1, . . . , i, . . . , P); then a trans- 
duction operation (Figure 9) is defined as follows: transduction 
is the transfer of genetic material from a Petri dish and bac- 
terium donors (p D , b D ) to a Petri dish and bacterium recipients 
(p R , b R ). When the transference involves a chromosome seg- 
ment, the result is a recombinant chromosome in the recipient 
Petri dish p R . However, the transference of a complete chro- 
mosome results in the substitution of one chromosome of the 
recipient Petri dish p R with the transferred one. It is impor- 
tant to note that "bacterium" and "Petri dish" terms are used 
throughout the paper as "chromosome" and "sub-population" 
synonyms, respectively. Transduction requires the selection of 
the Petri dish and bacterium donors (p D , b D ), as well as the 
Petri dish and bacterium recipients (p R , b R ). In the reference 
(Perales-Gravan et al., 2013) we describe how transduction was 
conducted. 

The current PETRI algorithm (Figure 9B) uses a population 
size of AT, performing r e replicates, with P being the total number 




{p D ,t> D ) 



r 




(p D ,t> D ) 



FIGURE 9 | Transduction experiment. The figure shows transduction from 
donor Petri dish (p D ) and bacterium (b D ) to recipient Petri dish (p R ) and 
bacterium {b R ). In the figure, P is the total number of Petri dishes (or 
sub-populations), N is the number of bacteria (or population size) per Petri 
dish and r e the number of experimental replicates. 
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of Petri dishes or sub-populations. Thus, we performed a num- 
ber of r e . P trials of each simulation experiment. The algorithm 
cycles through epochs, searching for an optimum solution until 
a maximum of G generations is reached. Once (p D , b D ) and 
(p R , b R ) are selected, only one "bacteriophage" is assumed to 
participate during each transduction event. The PETRI algorithm 
is summarized in the following pseudocode description: 

/* PETRI: Genetic Algorithm with Transduction */ 

1. t:=0; 

2. Initialization: Generate P Petri dishes (or sub-populations) 
with N random bacteria (or chromosomes). 

3. WHILE not stop condition DO 
/* Genetic Algorithm */ 

(3.1) FOR each P Petri dish DO 
Evaluation of chromosomes 
Selection 

Conjugation or Crossover (one-point, two-point) 
Mutation 
(3.2) END FOR 

/* End of Genetic Algorithm */ 

4. Transduction: (pD, bD) (pR, bR) 

5. t:=t+l; 

6. END WHILE; 

/* End of PETRI */ 

We studied the performance of the simulated transduction by 
considering three optimization problems that are described in 
sufficient detail in Perales-Gravan et al. (2013). The first problem 
uses a benchmark function, the second one is the 0/1 knapsack 
problem, and finally we illustrated the usefulness of transduction 
in the problem of designing an AM radio receiver (Perales-Gravan 
and Lahoz-Beltra, 2008). 

CONCLUSIONS 

In this chapter we have reviewed some of the models, simula- 
tions and theories that we have been working in recent years. 
The main conclusion of our work is that the bacterial cell can 
be seen as a form of natural computing (Lahoz-Beltra, 2012b), 
to which we have referred to in this chapter as "bacterial comput- 
ing." In bacteria computing capability emerges from two related 
processes: learning and evolution, being illustrated in this chap- 
ter several examples of hardware inspired in these processes. The 
possible impact of bacterial computing is not only to show how 
evolvable hardware can be used as a modeling framework in the 
simulation of learning and evolution in bacteria, also it promotes 
how electronic circuits could be designed based on "bacterial 
algorithms." 



In future, there will be wide range of applications of bacte- 
rial computing. For instance, recently Ran et al. (2012) designed 
a device made of DNA inserted into bacterial cell that works like 
a diagnostic computer. This molecular device works like a NOR 
logical gate being programmed to check for the presence of two 
transcription factors in such a way that responds by creating a 
protein that emits a green visible light — a sign of a positive diag- 
nosis. Also, a few year ago a new technique has been developed 
to save data in bacteria. According to Yachie et al. (2007) up to 
100 bits of data can be saved in each microorganism. Scientists 
successfully encoded and saved the phrase "e = mc 2 1905" the 
DNA of Bacillus subtilis. However, at present bacterial computing 
is a branch of synthetic biology. At present, bacterial comput- 
ing depends of synthetic biology and the latter is still in its early 
stages. In order to design and build genuine "bacterial computers" 
we will need a better understanding of the operation of complex 
biological systems, i.e., bacteria. 
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